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Abstract 

An effective two-body interaction is constructed from a new Reid- 
like NN potential for a large no-core space consisting of six major 
shells and is used to generate the shell-model properties for light nuclei 
from A=2 to 6. For practical reasons, the model space is partially 
truncated for A— 6. Binding energies and other physical observables 
are calculated and they compare favorably with experiment. 
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1 Introduction 



In traditional nuclear shell- model calculations, only a few particles or holes 
with respect to a closed shell are treated as active within a restricted model 
space. In a well-studied example, ^^O, the model space contains one ma- 
jor shell, the Is-Od shell, with two valence nucleons. These calculations 
require effective-interaction matrix elements along with calculated or empir- 
ical single-particle (s.p.) energies as input. The effective interaction could 
either be "phenomenological" (see, for example, Refs.|l], ^, ^]) or "realis- 
tic" (see, for example, Refs.[Q, ^]), depending on how it is obtained. Both 
types of effective interactions have been substantially used with success, 
when good agreement with experimental spectra is taken as the criterion. A 
phenomenological interaction might be obtained by fitting selected energy 
spectra and electromagnetic properties of the nuclei in the region of interest. 

In the case of a realistic interaction, which is our main concern in this 
work, one usually starts with the Brueckner reaction matrix G Q (i.e., lad- 
der diagrams) calculated from a realistic nucleon-nucleon (NN) potential 
for a specified model space, and evaluates other diagrams (e.g. folded di- 
agrams and/or core-polarization diagrams Q) using the G matrix. This 
renormalization procedure is incomplete, however, because so far, the core- 
polarization diagrams can only be included to, at most, the third order in the 
perturbation-theory expansion |9|, [l^. The incompleteness here presents a 
serious problem because convergence has not been attained within the lowest 
few orders of the perturbation expansion 11|. Similar uncertainties exist 
when calculating the effective operators to be used in the model space. 



Recently, attempts |13|, |15|, |16|, |1^] have been made to avoid the 
above difficulty by adopting a no-core model space in which all nucleons 
in a nucleus are treated as active. It is considerably simpler to obtain the 
effective interaction for a no-core space since there are no hole lines and 
the complicated core-polarization processes are absent. Consequently, one 
is left with the ladder diagrams and the folded diagrams for the effective 
interaction which may include effective many-body contributions. 

Within the concept of no-core calculations, it is important to distinguish 
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two cases. In the case we call a "full" no-core calculation, one selects a set 
of d model-space s.p. states and then generates the configurations where 
all nucleons can occupy all orbitals in all possible, Pauli-allowed, ways. In 
an "iVmax^^^ truncated" no-core calculation, only those configurations are 
retained from the full no-core case in which there are up-to and including 
-^max^^ excitations of the lowest unperturbed configuration (in harmonic 
oscillator notation) of the A nucleons. 

To be more specific, let us consider a 2hn (i.e., iVinax=2) shell-model 
diagonalization for ^Li in a no-core model space consisting of the lowest four 
major shells (Os, Op, Is-Od and Ip-Of). In this case, the configuration with a 
hole in the Osi shell and a particle in the Is-Od shell [i.e. (Os)^(Op)^(lsOd)^] 
is allowed. The configuration [i.e. (Os)^(Op)^] is also taken into account. 
However, one cannot claim to have performed a "full" no-corc calculation 
because only one or two, out of four Os nucleons, are allowed to be excited to 
the higher shells. Namely, in this 2hO, truncated calculation, not all nucleons 
are active, and there still is a partially inert core. 

On the other hand, if one includes 2s-ld-0g and 2p-lf-0h major shells 
and performs a Afiil calculation for the same nucleus, the configuration (Op)^ 
will be allowed, leaving no nucleons in the "core" orbital Osi. Although 

2 

such a calculation is still restricted, it is surely more complete than the 2hO, 
calculation. It is not currently possible to actually carry out a full no-core 
calculation in many cases we want to investigate. Our hope is that as A'^max 
increases, the results will converge and approach those of the full no-core 
calculation. 

Another practical issue of working with an N^s^^^O, truncated no-core 
calculation is that this facilitates the accurate treatment of the spurious 
center-of-mass (cm.) motion. If (Nq + l)hO, is defined as the minimum s.p. 
excitation energy needed to lift a nucleon to the lowest unperturbed state 
outside the model space (A''o=4 for the ®Li example above in the model space 
through the 2p-l/-0/i shell), then, for A^max = Nq, it is possible to obtain 
no-core shell-model wavefunctions free of spurious cm. motion. 

It is an ultimate goal of the nuclear shell model to be able to start with 
a realistic NN potential and obtain unambiguous and converged results 
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against the changes in the size of the model space and in the choice of the 
unperturbed basis. Convergence with model-space size means convergence 
with increasing A^n^ax and increasing d — a dual convergence criterion. 

Encouraging results have been obtained recently for very light nuclei in 
Ref . , by Ceuleneer et al. They have performed a shell-model calculation 
for the T=0 states in ^He, where up to IQfiQ. excitations from the basic 
configuration [(Osi)^] are allowed. The only input to the calculation is a set 
of two-body matrix elements (tbme) of a modified Sussex interaction. Since 
this effective interaction does not have a theoretically derived model-space 
dependence, they multiplied all two-body matrix elements by a model-space 
dependent parameter which is adjusted to get the correct binding energy. 
The step of deriving the dependence of this parameter on the model space 
size is now required to complete the dual convergence test. 

In this work, we will adopt a large no-core harmonic-oscillator (HO) 
model space consisting of 6 major shells (from Os to 2p-lf-0h). We will 
consider several light nuclei ranging from A=2 to A=6. An effective inter- 
action will be constructed for the above model space from a new Reid-like 



A^A^ potential (Reid93) provided by the Nijmegen group |19]. Note that we 
will use effective interactions constructed in exactly the same manner for 
all the nuclei considered here. We will follow an approach that favors the 
more accurate treatment of the spurious cm. motion and attempts to min- 
imize the neglect to the two-body "ladder" scattering procedures. We have 
designed an even more accurate approach along these lines which involves 
excitation-dependent effective interactions and will be reported in a future 
work 

It is an established fact that for a small model space, a mass-dependent 
two-body effective interaction gives an overall better description ||2|, But 
we anticipate that such a mass dependence will become weaker as the size 
of the no-core model space is increased. Similarly, we expect the effective 
many-body forces to decrease with a increasing model space. Indeed, if an 
infinitely large model space is used, the effective interaction reverts back to 
the A^A^ potential v, whose matrix elements are clearly independent of the 
nucleus under consideration. Throughout the remainder of this work, we 
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assume the model space is sufficiently large and the s.p. basis is sufficiently 
realistic that we can neglect the effective many-body interactions. This will 
be investigated in a future effort, which also addresses the rate our methods 
approach the goal of satisfying the dual convergence criteria. 

2 Effective Interaction 

For a no-core model space, the core-polarization diagrams are not present, 
and the two-body effective interaction is simply the G matrix plus the 
folded diagrams series [^. The G matrix is the sum of the ladder diagram 
series which represents the multiple scattering processes of two nucleons in a 
nuclear medium. We continue to follow our philosophy given in Ref.[|l^] for 
the no-core G-matrix in large spaces which treats two-particle scattering via 
a realistic A^A^ interaction vi2 in an "external" field, u, which is provided by 
the other nucleons in the same nucleus. Thus, we write 



where h = t + u is the one-body Hamiltonian and u is the nuclear mean field. 
The quantity uj is the starting energy, which represents the initial energy of 
the two in-medium nucleons. The Pauli operator Q excludes the scattering 
of the two nucleons into the intermediate states which are inside the model 
space. It is therefore related to the choice of the model space and will be 
specified in the next section. 

A rigid prescription for the nuclear mean field u is not necessary since 
the results will be independent of u once the dual convergence criteria are 
satisfied. In most practical calculations, u is approximated by a one-body 
potential of a simple and convenient form. The two most common choices 
for u are a shifted HO potential and zero: 




Vl2 + Vi2 



Q 



(1) 





(2) 
(3) 



u{r) = 0. 
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The latter choice corresponds to a plane- wave basis. Some hybrid approaches 
have been developed which use oscillator states for low-lying orbitals and 
plane waves, orthogonalized to the oscillator states, for all the remaining 
states l2|, |4|. 

Although a shifted HO potential (||) does not have the expected asymp- 
totic behavior of vanishing exponentially at large r, it was argued in Ref.|17] 
that the shape of the assumed u at large r might not be very important 
since, except for some weakly bound states, nucleons are unlikely to move 
far beyond the nuclear mass radius. 



One may further notice from Refs.[Q, 25 1 that the two seemingly very 
different one-body potentials in Eqs.(|2|,^) actually led to rather similar G- 
matrix elements, provided one makes a careful choice for the starting energy 
(related to the choice of u) . Note that the constant shift Vq in Eq. (^) is more 
a matter of convenience, as a shift of 2Vq can be made in the starting energy 
uj in Eq.(|l]), i.e., 

uj = uj'-2Vo, (4) 

to cancel out Vq in the energy denominator [[T^ ]. 

In this work, we will approximate the nuclear mean field by the HO po- 



tential not only because this simplifies the G-matrix calculation |25, 2£], 
but, more importantly, for the reason that this makes possible an exact 
removal of the effects of the spurious cm. motion from our many-body 
wavefunctions. Once the G matrix G{lo) is obtained as a function of the 
starting energy, it will not be difficult to evaluate the folded diagrams us- 



ing the techniques developed by Kuo and Krenciglowa [27| and by Lee and 



Suzuki 1 28 1 and to obtain a starting-energy-independent effective two-body 
interaction (denoted by v^g). 

(2) 

One must bear in mind that v^g obtained in this procedure depends on 
the assumption made for the one-body potential in the G-matrix calculation. 
Especially in cases when the model spaces are small and we are further 
from satisfying the dual convergence criteria, it is important to use a u that 
sensibly represents the nuclear mean field so as to minimize the neglected 



many-body interactions [29| and higher than linear order "-ti" insertions. 
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In Ref.|17] it is shown that f^jj can be well approximated by the G matrix 
calculated at starting energies which are related to the initial unperturbed 
energy of the two nucleons in the ladder scattering processes in a simple way: 

= w + 2Vb = ea + efe + A, (5) 

where e = (2n^. + I + 3/2)?iil are the HO s.p. energies (a and b are the s.p. 
states that the two nucleons initially occupy). Such a state-dependent choice 
for Lo' will lead to a non-hermitian G matrix, but the non-hermicity is found 
to be small. The quantity A signifies the interaction energy between the two 



nucleons. In a specific application to Li, it has been shown |17] that for 



?iO=18 MeV, a value of -21 MeV for A results in G{u>') which is an excellent 
approximation to . 

In this work, we will follow Ref. ||l7[| and adopt the average of G{uj') and 
its conjugate calculated for a HO basis with ?iil=14 MeV with uj' given by 
Eq.@ as our approximation to . The parameter A is chosen to yield 
the experimental binding energy. Initially, one might expect that different 
values of A have to be used for different nuclei. But, quite surprisingly, we 
find that good agreement with experimental observables can be obtained 
with a nucleus-independent value of A (-35 MeV). 

Our shell-model Hamiltonian will now be written as 



/A \ A 

HsM = I X! ~ Tcm. I + X! + ^Coulomb + 
\i=l / i<i 



hn), (6) 



where U = -pf / (2m) are the one-body kinetic energies, Tc.m. = (Z^i Pi)^/ (2m^) 
the cm. kinetic energy and Vcouiomb the Coulomb interaction. The proton 
and neutron masses are taken to be the same. The last term (with A=10) 
in the above equation forces the cm. motion of the low-lying states in the 
calculated spectrum to be in its lowest HO configuration. 

We remark that our calculations involve no free parameters other than 
those used in calculating the G matrix, fiVt and A. Moreover, these two 
parameters are fixed at 14 MeV and -35 MeV, respectively, for all nuclei 
considered in the present work. 

We emphasize that in a no-core calculation, we are attempting to derive 
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all shell-model properties from an underlying Hamiltonian, H^m- Thus, 
there are no phenomenological s.p. energy terms in H^m- 



3 Results and Discussions 

As previously mentioned, we use a no-core model space containing the lowest 
six HO major shells with 7iO=14 MeV. For A < 4, we allow all OhO, through 
Ihfl configurations within the model space. For A > 4, we allow all OhO, 
through 5hi} configurations. Therefore, different Q operators have to be 
used in Eq.([l|) for A < 4 and A > 4: 

For A < 4 : Q = 1 for ni > 6, ris > 6, or ni + ?i2 > 8, (7) 

= otherwise; 

For A > 4 : Q = 1 for ui + > 6, (8) 

= otherwise. 



In the above equations, n=2nr+l is the principal quantum number for the 
HO s.p. states. It starts from with n=0 representing the first major 
shell (Os). For A=6, due to the computer memory limitation, the n=5 shell 
contains only the p orbitals; the / and h orbitals are left outside the model 
space. 

The shell-model matrix diagonalizations are performed for the Hamil- 
tonian HsM in Eq.(^) using the Many-Fermion-Dynamics code [ pO| . The 
results are given in Table I, which we will discuss in the following subsec- 
tions. The experimental results given in Table I are taken from Ref . [^] for 



A=3, Ref.g for A=4 and Ref.|33| for A=5 and 6. 



3.1 Binding Energies 

It is possible to obtain exact or nearly exact results for ground-state energies 



of the lightest nuclei by solving the Schrodinger or Faddeev |34| equations 
for realistic NN interactions. This has been done for the ground states of 
^H, ^He and ^He [35, 36, p^, pq]. Even for ^He and ^Li, preliminary results 



obtained by Wiringa in variational Monte Carlo calculations have appeared 
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pj9|| . Unlike the few-body approaches in which one obtains almost exact 
results, at least for the ground state, the effective-interaction shell-model 
approach involves some uncertainties due to the truncation of the space 
and the approximation made in calculating the effective interaction for the 
truncated space. Consequently, the shell-model approach to the above light 
nuclei, although being able to calculate for a given set of quantum numbers 
the excited states as easily as the lowest state, cannot match the few-body 
approach in the accuracy of the results for the ground state. Nevertheless, 
our ultimate goal is to satisfy the dual covergence criteria. By using a large 
no-core model space along with a reasonable effective interaction, we hope to 
demonstrate that, in spite of its present limitations, the effective Hamiltonian 
approach gives a useful description of the low-lying states in light nuclei. 

Our results are encouraging as can be seen from Table I. The calcu- 
lated binding energies for the deuteron, triton, ^He, ^He and ^Li are 2.103, 
8.589, 28.757, 25.960 and 30.648 MeV, respectively, agreeing quite well with 
the corresponding experimental values of 2.225, 8.482, 28.296, 27.410 and 
31.996 MeV. Of course, it is more relevant to compare our results to those 
obtained in the exact few-body approaches using the same potential. These 
approaches show that existing realistic A^A^ potentials underbind light nu- 
clei with A > 2. Our calculations involve a free parameter A which has 
been fixed at -35 MeV for all the nuclei considered. For A > 2, we can ob- 
tain smaller binding energies (in better agreement with exact calculations) 
by decreasing A (i.e., making it more negative) since the binding energies 
decrease monotonically with the decreasing A. Our adoption of a A value 
that fits experimental binding energies stems from an assumption that our 
neglected effective many-body forces and other corrections can largely cancel 
the neglected (and largely unknown) true many-body forces. 

It is worth mentioning that for the two-body system, the deuteron, it is 
possible to obtain exact results even with the effective-interaction approach 



1 41]. Our present results for the deuteron are not exact due to our ne- 
glect of the processes which are higher order in u. Our effective interaction 
nevertheless gives a reasonable deuteron binding energy. In section 3.3, we 



will further show that the deuteron magnetic dipole and electric quadrupole 
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moments also come out well. 



3.2 Excitation Spectra 

For the deuteron and the triton, we obtain only one bound state in the 
calculations, agreeing with experiment and with exact calculations. For the 
deuteron, the lowest state in the continuum is a J'^=0~^, T=l state, which 
is unbound by 1.65 MeV (i.e. 3.75 MeV above the ground state). For the 
triton, the lowest excited state is a J^=^ , T=^ state, unbound by 4.13 
MeV. It has a nearly degenerate J'^=^ , state, unbound by 4.28 MeV. 
The 7"= I states are even higher in energy. Therefore, these results do not 
support nn, pp, nnn or ppp bound states. 

For ^He, the experimental level sequence of the low-lying negative-parity 
states is correctly reproduced. The excitation energies are consistently higher 
than the experimental results [^] by about 2 to 3 MeV. These results are 
clearly an improvement over those obtained in our previous study [|l6|. In 
that study, the excitation energies of these same states were obtained in a 
smaller model space, including only four major shells, and were found to be 
as much as 6 MeV too high when compared with experiment (see Table I in 
Ref . [|l^] ) . The better results we obtain here should be attributed mainly to 
the larger model space and the improved NN interaction. From a theoretical 
viewpoint, we have also improved the G matrix by using a state-dependent 
starting energy of Eq.(^) (rather than at a constant starting energy as in 



Ref. [p!q] ) which better approximates the full effective interaction |1?]. 

We obtain the first excited state (J'^=0+, T=0) in '^He at an excitation 
energy of 26.135 MeV. This is about 6 MeV higher than experiment but it 
is about 7.7 MeV lower than the previous result (33.807 MeV) for the four- 
major-shell space [^B]. Again, the larger model space used in this work is 
largely responsible for the decrease in energy. A more accurate description 
of this state will require an even larger space. Indeed, in Ref. |T^ where a 
modified Sussex interaction is used, excellent agreement with experiment is 
obtained for this state only when up to lOhil. configurations are included. 

The calculated ground state in ^He is dominated by the (Os)^ configu- 
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ration but it has a considerable amount of "Ip-l/i" configuration (Os)^(ls)^ 
and "2p-2/i" configuration (Os)^(Op)^. The state is dominated by the 
(Os)^(ls)^ configuration while the (Os)^(Op)^ and (Os)^ components are also 
quite significant. 

The occupancies of the various model-space orbitals in the 0^ and 0^ 
states are 

|4He:0+) - (Os)3-^25(Op)0.i60(^^)0.200(Q^j^g^gj^gllg)0.ii4. (9) 
|^He:0+) = {Osf '''^{Opf ''^\sd)°-'^^ (other shellsf-'''^ (10) 

Relative to the ground state, the O2 state has only about 50% of the "breath- 
ing mode" (Os)^^(ls)^. However, this result depends on the choice of the 
s.p. basis. If, for example, a Hartree-Fock basis were used, the oscillator- 
basis s-states would mix to produce the HF s-states (e.g. Oshf and Ishf 
etc.) so that the admixture of the (Oshf)^(I'Shf)^ component in the ground 
state would likely be much smaller. This would lead to a larger amount of 
(Oshf)~"'^(1'Shf)"'^ in the 0^ state. 

Note that although our model space is not sufficiently large to reproduce 
the first excited 0"*" state in ^He at the experimental energy, it does a fairly 
good job for the "ITiO" states. This gives us confidence in the results for the 
low-lying states in ^He, which we present below. 

The first excited state in ^He {J'^=-^ , T=^) is obtained at an energy of 
3.112 MeV, within the range of 3 to 5 MeV given in Ref . p3[| . The low- lying 
positive-parity states are also of interest. Experimentally, there is a famous 
J'^=^~^ state at 16.75 MeV. It has a dominant (Os)'^(Op)^ configuration and 
can be thought of as the ground state of ^Li coupled to a Os hole. This state 
corresponds to our calculated 20.445 MeV state in Table I, which has the 
following occupation probabilities: 

|20.445MeV : = {Osf-^^^iOp)^-^°'^{sdf-^^'^{othei shells)°-°9\ (11) 

22 

to be compared to the occupation probabilities for the ground state of ®Li 
I^Li : 1+) = (0s)3-63i(0p)2-06i(^rf)0.222(^^i^g^ shells)"-^^^ (12) 

Note that the fact that the calculated energy of this J^=^ state is about 3.7 
MeV higher than experiment is more or less consistent with what we have 



11 



seen in the case of ^He where the "IhQ" low-lying negative-parity states 
came out about 2 to 3 MeV higher than experiment. It will be interesting 
to track the energies of these states as well as the excited states in ^He with 
increasing model-space size. 

Above the two s.p. states (the ground state and the first excited ^-^ 
state) and below the 20.445 MeV |^ state, our calculation also gives three 
positive-parity "ITiO" states, a ^ state at 7.437 MeV and nearly degenerate 
§ and I states at 14.206 and 14.439 MeV, respectively. These states are 
dominated by the configurations (Os)^(sci)^ and (Os)^(Op)^. The occupancies 
of the orbitals in the model space are: 

1 + 

|7.437MeV:-) = (Osf'^^^(0pf-^^^(sd)°-^^\othershelh)°-^^'^; (13) 
2 1 

|14.206MeV: - (0s)3 i97(0p)°-^28(^^)0.67i(o^j^g^ gj^gllg)0.204. (^4) 

2 1 

I14.439MeV: = (Osf -^'^^iOpy ^^^sdf -^^^ {other shelh)°-^^^ . (15) 

2 1 



There have been previous theoretical predictions |33, 42] that there is a 
state at about 5 MeV and |^ and |^ states at about 12 MeV. These 

predictions have not been fully confirmed experimentally, but they are well 

supported by our results, again, keeping in mind that our calculated "l?ir2" 

states are probably about two or three MeV too high. 

In addition to the above low-lying states, we have also listed in Table I a 

few other bound states of ^He which have an energy not much higher than 

the experimental 16.75 MeV state. 

The low-lying energy spectrum of ^Li obtained in this calculation does 



not show much improvement over that in Ref.|16|. It again appears to be 
more spread-out than the experimental spectrum. 



3.3 Ml and E2 Moments 

Since we are using a large no-core model space, we choose to use bare op- 
erators (ep=l, 6^=0, (7p=5.586, (7^ =-3. 826, gp=1.0, §1^=0.0) to calculate the 
magnetic dipole (Ml) and electric quadrupole (E2) moments in leading ap- 
proximation. The calculated results are also given in Table I. It should be 
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emphasized that only the nucleonic degrees of freedom are taken into account 
in calculating these moments. Proper considerations have to be given to the 
effects of the meson exchange currents (MEC) before critical conclusions can 
be drawn from the comparison of the calculated moments (especially the Ml 
moment) with data. 

The calculated Ml moment fi for the deuteron is 0.857^Ar. This agrees 
with the experimental result of 0.8574/XAr. However, this fortuitous agree- 
ment will be vitiated to the extent that the ignored MEC contribution is 
significant. Even if the MEC effect is negligible, the value that we obtained 
for the deuteron Ml moment is not theoretically exact. This is made evident 
in the discussion below. 

The deuteron Ml moment is related to the I?-state probability Pd as: 

/i(2H) = Ps^^i^Si) + PdK^Di) = (1 - Pd)0.880 + PdO.310 (/x^). (16) 
With this equation, a calculated value of 0.857/iAr for //(^H) leads to Pd=4.0%. 



However, the exact Pd for the Reid93 potential is in fact 5.7% [40|, imply- 
mg a /i(2H) of 0.848;Uiv. We, therefore, see that the tensor force is somehow 
weakened when we go from the bare A^A^ potential to the effective shell- 
model interaction in Eq.(^) for our no-core model space. This infers the 
size of the neglected contribution to the magnetic moment operator arising 



in the theory of effective operators. It has been shown in Ref.|43] that the 
tensor force strength can be further reduced by core-polarization diagrams 
(mainly the Bertsch bubble diagram [Q]) that one must take into account 
when calculating the effective interaction for a small, one-major-shell, model 
space outside an inert core. 

The calculated deuteron quadrupole moment Q is 0.242efm^, somewhat 
smaller than the experimental value of 0.286efm^. This agrees with the 
above observation that the effective tensor force in our no-core shell-model 
interaction is weaker than that in the original NN potential. The reduced 
quadrupole moment may also arise from the fact that its operator involves 
a radial dependence (r^) which needs to be renormalized when we truncate 
the infinite Hilbert space to our finite-size no-core HO model space. Thus 
we reason that, for our model space, the renormalization effects are larger 
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for the E2 operator than for the Ml operator which does not have a radial 
dependence. 

The need for using an effective operator to evaluate the root-mean- 
squared (rms) radius (or any other observable that involves it) is evident 
from Table I, where the calculated rms point radius \J {I'p) for the proton in 
the deuteron is 1.488 fm, significantly smaller than the experimental value 
of 1.95 fm. The large renormalization of the rms radius operator required 
for the deuteron is not surprising since it is a very loosely bound system, the 
wave function obtained in the truncated HO model space does not represent 
the exact wave function very well. The calculated \J {fp) value for ^Li is also 
smaller than the experimnetal value. However, the results of \J {r^) for '^H 
and ^He are in good agreement with experiment. Note that we have evalu- 
ated these rms radii with "intrinsic" wavefunctions so the quoted results are 
free of spurious cm. contributions. 

Our calculated Ml moment for the triton is 2.659/iAr, about 11% smaller 
than the experimental value of 2.979//iv- To a large extent, this discrepancy 
may be explained by the MEG effects that we have not taken into account. 
Indeed, in Ref.[^], it is shown that the inclusion of the MEG effects in a 
model-dependent way leads to a 14%'s increase in the triton Ml moment 
from 2.588^Ar to 3.010/iAr, in close agreement with experiment. 

For the ground state of ^He, the calculated Ml and E2 moments are 
-1.864^7V and -0.332efm^, respectively. Again, the MEG effects have to be 
considered when comparing these results with experimental data, which, to 
our knowledge, are not available. 

It has been difficult in the past for theory to reproduce the E2 moment for 
the ground state of ^Li. However, the calculated E2 moment is -O.llGefm^, 
which is remarkably close to the experimental value of -0.082efm^. Our 
calculated Ml moment is O.SSI^at, which is about 3.5% higher than the 
experimental result of 0.822/ijv. 
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3.4 Effects of the Coulomb Interaction 

Since we include the Coulomb interaction, the isospin symmetry is not 
strictly conserved. But the isospin impurity caused by the Coulomb in- 
teraction is generally very small. For the bound states of '^H and ^He, the 
calculated values for isospin 



tcalc = ^ (17) 

are 0.500000 and 0.500022, respectively. Note that has only one proton so 
isospin is still a good quantum number. In ^He, the calculated isospin shows 
only a 0.0044% deviation from the half- integer value. T^aic is 0.000046 for 
the ground state of ^He; it is 0.500016 and 0.500024 for the ground states of 
^He and ^Li, respectively. The small isospin impurity for the ground states 
in these nuclei is due to the fact that all these states do not have any nearby 
state with the same J'^ but a different T. From perturbation theory, one 
knows that the relatively weak Coulomb interaction will not induce much 
isospin mixing to these isolated states. 

The Coulomb interaction has sizable effects on the absolute energies of 
the system, as is well known. Our calculation shows that, due to the Coulomb 
repulsion, the binding energy of ^He is 0.725 MeV less than that of and 
the binding energy of ^Li is 1.024 MeV less than that of ^He. The experi- 
mentally observed differences in the binding energies for the above two pairs 
are 0.764 and 1.073 MeV, respectively. They are quite close to our calculated 
values, as one might expect since the Coulomb interaction is a perturbation 
in these light systems. Nevertheless, our results for the Coulomb energy are 
model-dependent [in that the matrix elements of the Coulomb interaction 
in the shell-model Hamiltonian were evaluated using a HO basis and 
possible renormalization corrections from the excluded space were ignored]. 
A smaller Coulomb effect of about 0.74 MeV was obtained in a more model- 
independent analysis |45| for the ^H-^He pair. It is also believed that other 



charge-symmetry breaking effects contribute to the difference between the 
binding energies of and ^He as well 
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4 Conclusions 



In this work, we have constructed an effective interaction for a six-major-sheh 
no-core model space from a new, Reid-Hke, NN potential (Reid93) from the 



Nijmegen group |19]. The effective interaction has been applied to calculate 
nuclear structure properties for a few light nuclei, ranging from the deuteron 
to ^Li. The results are very encouraging. Not only are the binding energies 
of these nuclei well reproduced, the energy spectra are also in good agree- 
ment with experiment. In particular, the experimental level sequence of the 
low-lying negative-parity states in ^He is correctly reproduced, although the 
excitation energies are about 2 to 3 MeV higher than experiment. Based on 
our current and previous efforts, we expect that this discrepancy will be re- 
duced as we more closely satisfy the dual convergence criteria — convergence 
against increasing A^'max and d, where A^'max signifies the highest unperturbed 
energy of the configurations taken into account and d represents the number 
of s.p. states included in the model space. 

The magnetic dipole and electric quadrupole moments, calculated using 
bare operators with meson-exchange-currents effects neglected, are also in 
reasonable agreement with experiment. 

For ^He, in addition to the two low-lying s.p. negative-parity states | 
and ^ , we have obtained a low-lying state at about 7.4 MeV and two 
nearly degenerated states (|^ and |^) at 14.2 MeV and 14.4 MeV. The latter 
three, dominated by the configuration (Os)'*(s(i)^, are mainly s.p. states with 
one {sd) neutron coupled to the ground state of ^He. The actual energies 
of these predominantly "l?iil" states could be about a few MeV lower, as in 
the case of ^He. The previous theoretical predictions of a state at about 
5 MeV and | and | states at about 12 MeV are therefore well supported 
by our results. The 16.75 MeV state, resulted from the ground state of ^Li 
with a (Os) proton removed, is reproduced at an energy of 20.445 MeV. 

The Coulomb interaction, which is included in the calculations, accounts 
for the bulk part of the differences in the experimental binding energies of 
mirror pairs (^H-^He and ^He-^Li). We have also seen that the Coulomb 
interaction induces a very small amount of isospin impurity to the ground 
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states of the light nuclei considered. 

An extension of the current approach to heavier Op-shell nuclei will be 
straightforward. Our results for A=2 to 6 have given us optimism that our 
approach would be able to give a good description of neighboring nuclei as 
well. This is presently being investigated. 

Of course, since the size of the shell-model matrix increases quite dra- 
matically with the increasing number of nucleons, it is unlikely at the present 
time that one can apply the no-core approach to a much heavier nucleus, like 
^'^Ca. In this regard, the Monte Carlo shell-model approach |47], in which 
the size of the calculations increases only moderately with the number of 
active nucleons, offers some promise. 
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Table I. The results for ^H, '^H, ^He, ^He and ^Li obtained in large no-core (con- 
sisting of 6 HO major shells) shell-model calculations. The experimental data are 
taken from Refs.|31, p2l B3|. In the Table, Eb is the binding energy (in MeV); 



Ex{Jn,T) the excitation energy (in MeV) of the J^jT state. The ground-state rms 
point radius for protons ^ (r^) (in fm), electric quadrupole moment Q (in efm^) and 

magnetic dipole moment fi (in /ijv) are also listed. The "experimental" ^ (r^) was 
deduced from the charge rms radius \/Jr^ through (ignoring the neutron charge 
distribution and other higher-order effects and assuming a proton rms charge radius 
of 0.81 fm) (rl) = (r^) - O.Sl^. 
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") Low-lying positive-parity states (e.g. a ^'^^^ , T=^ state at ~5 MeV and J^—^ , 
r=i and J^=|^, T=i states at '-1 2 MeV) are predicted to exist. See Ref.[||| for 
more details. 

*) We identify the calculated 20.445 MeV state as the experimental 16.75 MeV state, 
because the calculated state is dominated by the (Os)^(Op)^ configuration. 
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